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Description 

[0001] The present invention relates to a method according to the preamble of claim 1 and to an apparatus for 
implementing said method. 

s [0002] Such a method is known from WO-A-9000367 in which a system for the management of the hemodynamic 
state of a patient is disclosed. With this respect electrocardiogram signal data from a patient and in particular a time- 
derivative thoracic impedance signal is processed. 

Background of the Invention 

10 

[0003] Cardiac monitoring to which the present invention relates includes both the determination of heart rate (HR) 
from electrocardiogram (EKG) signals and the determination of heart stroke volume (SV) from thoracic impedance 
signals, from which cardiac output (CO) can be estimated. 

[0004] Heart rate can be determined in a number of ways. The phonocardiogram is considered among the most 
15 accurate methods of determining heart rate. However, due to the practical difficulties in using it, the phonocardiogram 
method is generally not employed for any continuous or long-term monitoring. 

[0005] Heart rate is most typically determined by the electrocardiogram (EKG). The analog EKG signal typically 
displays electrocardio events as perturbations referred to as waves. The heartbeat is most clearly reflected in the EKG 
signal as an R wave peak between a pair of adjoining Q and S wave valleys. The basic and commonly used method 

20 of automatically identifying the QRS wave pulses in point is the threshold method in which the rate of voltage change 
between consecutive data points of the EKG signal is monitored and compared with a threshold value. Slopes exceed- 
ing the threshold value are deemed to be associated with a portion of the QRS pulse. While this method commonly 
detects the interval between consecutive R waves successfully more than eighty percent of the time, it typically has 
difficulty in dealing with sources of irregular signal components such as pacemakers, muscle noise, 60 Hz interference 

25 as well as nearby T or P waves which may also be associated with significant slope changes. 

[0006] Hemodynamic monitoring of the heart can provide very valuable physiological information regarding the func- 
tional state of the myocardium, which is intimately related to its mechanical behavior. The quantitative measurement 
of blood flow, or the cardiac output (CO), is one of the most useful parameters in assessing cardiac capability, but it is 
also one of the most difficult to measure. It cannot be accomplished with the electrocardiogram (EKG) which does not 

30 reflect the real pumping action of the heart. 

[0007] Both invasive and non-invasive methods are available for measurement of cardiac output. The invasive meth- 
ods are considered the most accurate. The risks associated with them are often an unacceptable trade-off, for they 
require direct access to the arterial circulation. In addition, invasive methods are not suitable for repetitive measure- 
ments and normally cannot be performed outside a hospital. Furthermore, invasive methods are very demanding in 

35 terms of time consumption and the need for skilled personnel. 

[0008] Impedance cardiography has been found to be one non-invasive method with the potential for monitoring the 
mechanical activity of the heart with virtually no risk. It can be conveniently handled by nursing and non-technical staff. 
It can usually tolerate moderate patient movement and can be left unattended for continuous and long-term monitoring. 
[0009] U.S. Patent No. 3,340,867, now RE 30,101, to Kubicek et al. discloses an impedance plethysmograph which 

40 employs four electrodes, two around the neck and two around the torso of a patient, to provide an impedance difference 
signal from the two center electrodes. The outermost pair of electrodes apply a small magnitude, high frequency al- 
ternating current to the patient while the inner pair of electrodes were used to sense voltage levels on the patient above 
and below the patient's heart. The impedances of the patient at each of the inner pair of electrodes could be determined 
from the sensed voltages and known applied currents. 

45 [0010] According to Kubicek et al., stroke volume (SV) is related to impedance Z as follows: 

SV=R(LVZ 0 ) 2 (VET) (dZ/dt m[n ) 

50 where R is blood resistivity, L is the distance between the inner voltage sensing electrodes, 2^ is the mean thoracic 
impedance determined from the inner, voltage sensing electrodes, VET is the ventricular ejection time and dZ/dt^ is 
the maximum negative slope change of time-differentiated impedance signal, which is the time-differentiated difference 
between the impedances determined at the center two electrodes. The above equation is referred to as Kubicek's 
equation. Cardiac output per minute is stroke volume times heart rate in beats per minute. 

55 [0011] The Kubicek equation is based upon a parallel column model of the thorax in which it is assumed: 

(1 ) the thorax is a cylinder, consisting of two electrically conducting tissue paths, also of cylindrical form with uniform 
cross-sectional areas and homogenous conducting materials, the first path being the blood with a relatively low 
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resistivity and the second path being all other tissues with relatively high resistivities; 

(2) the relationship between the maximum impedance change and the stroke volume during the cardiac cycle is 
linear; 

(3) impedance measurements of the individual specific tissue volumes are not very useful in developing the model 
5 (the parallel columns model relies on the intact thoracic measurements) ; and 

(4) at 100 kHz frequency, a physiologically safe frequency, the relative thoracic impedance changes are at a max- 
imum, the effects of polarization are negligible, and the reactive component of impedance is small, especially when 
compared to the real component, so that the reactance could be ignored in determining impedance with only a 
small error. 

10 

[0012] Yet another model and equation for relating impedance and stroke volume has been proposed by Sramek. 
According to Sramek, stroke volume (SV) is related to impedance Z as follows: 

fff SV=[(0.17H) 3 /4;2]*[VET]*[dZ/dt min /Z o ] 

where H is the patient's height. The Sramek equation is based upon a frustoconical model of the thorax. The Sramek 
model illustrates some improvement and accuracy over the Kubicek model but its major assumptions are still similar 
to those of the Kubicek model. 

20 [0013] Despite its advantages, impedance cardiography has not been well accepted by clinicians for three primary 
reasons: 

(1) poor correlation of the methods and models with the results of the more accepted invasive techniques; 

(2) still a relatively high dependance on skilled technical operators; and 

25 (3) still a discomfort to and/or disturbance of patients associated with the use and application of band electrodes. 

It is believed that poor correlation, the primary reason, can be traced back to a single source, namely the continuing 
inability to relate impedance cardiography and its mathematical model directly to cardiac physiology. 
[0014] The following are definitions and abbreviations of some of the terms used frequently herein: 
30 [001 5] Heart Rate (HR): the number of times the heart contracts each minute. 

[001 6] Ventricular Ejection Time (VET): the time interval of the opening and closing of aortic value during which there 
is movement of blood out of a ventricle. 

[0017] Stroke Volume (SV): the volume of blood pumped out by a ventricle (in particular the left ventricle) with each 
contraction of the heart. 

35 [001 8] Cardiac Output (CO): the amount of btood pumped out of the heart into the systemic circulation each minute. 
[0019] Ejection Fraction (EF): the ratio SV/EDV, which is the percentage of blood in a ventricle ejected with each 
contraction; it is directly related to the strength of the heart with <50% considered abnormal. 
[0020] End Diastolic Volume (EDV): the volume of blood that fills the ventricle before ejection. 
[0021] It would be desirable to determine heart rate more accurately than can be determined using the cardiogram 

*o threshold method currently employed. 

[0022] It further would be desirable to provide non-invasive, cardiographic impedance monitoring to estimate stroke 
volumes, cardiac outputs and related cardiac function parameters which correlate more closely with the stroke volumes, 
cardiac outputs and the like determined by means of recognized, accepted invasive procedures, but which does not 
require of operators the technical skills required by current impedance cardiograph systems, and does minimize dis- 

45 comfort to the patient on which the system is used, thereby permitting relatively long-term monitoring of the patient's 
condition. 

Summary of the Invention 

so [0023] One aspect of the invention is a method of processing a time-derivative, thoracic impedance signal generated 
from a patient to identify events in the time-derivative impedance signal associated with beats of the patient's heart. 
The method comprises the steps of generating a time-frequency distribution of data from the time-derivative impedance 
signal for a selected portion of the time-derivative impedance signal spanning a single heartbeat of the patient and 
identifying a time of occurrence of at least one cardiac event in the selected portion of said time-derivative impedance 

55 signal from the time-frequency distribution. 

Preferred aspects of the invention are methods of processing electrocardiogram signal data from a patient or of deter- 
mining stroke volume of a patient's heart from thoracic impedance measurements. 

Another aspect of the invention is an apparatus for implementing the method of processing a time-derivative, thoracic 
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impedance signal generated from a patient The apparatus is characterized by the provision of a first pair of electrodes 
on opposing upper and lower sides of the patient's heart to which a fluctuating excitation current is applied, by the 
provision of a second pair of electrodes on opposing sides of the patient's heart between the first pair of electrodes in 
order to generate a first signal related to impedance of the patient at a first electrode of the second pair and a second 
5 signal related to impedance of the patient of a second electrode of the second pair, by means for identifying an imped- 
ance from each of the first and second signals occurring simultaneously with one heartbeat of the patient including 
means for determining the time-derivative of said signals and means for generating a time-frequency distribution of 
data from the time-derivative of said signals, and by means for multiplying together the identified impedances of each 
of the first and second signals to estimate a cardiac parameter of the patient's heart. 

10 

Brief Description of the Drawings 

[0024] In the drawings, like numerals are employed to indicate like elements throughout. 

[0025] The foregoing summary as well as the following detailed description of the preferred embodiments of the 
is invention will be better understood when read in conjunction with the appended drawings. For the purpose of illustrating 
the invention, there is shown on the drawings embodiments which are presently preferred. It should be understood, 
however, that the invention is not limited to the precise arrangements and instrumentality shown. 

Fig. 1 depicts the cardiography monitoring system of the present invention with the electrodes thereof applied to 
20 a patient P; 

Fig. 2a depicts a portion of an exemplary EKG signal spanning a single heartbeat; 

Fig. 2b depicts a portion of an exemplary inverted dZ/dt signal spanning a single heartbeat; 

Fig. 2c depicts a portion of an exemplary phonocardiogram signal spanning the same heartbeat shown in the d27 

dt signal; 

25 Figs. 3a-3b depict in block diagram form, the steps followed by the processor of the system in determining the 

various cardiac parameters identified and calculated by the system; 

Fig. 4 depicts schematically a linear phase delay filter for combined low-pass, high-pass filtering; 
Fig. 5 depicts in block diagram form the detailed steps of the EKG signal processing; 

Fig. 6 depicts in block diagram form the detailed steps of the time differential thoracic impedance signal processing; 
30 and 

Fig. 7 depicts a two-dimensional representation of local frequency change determined from an exemplary Spec- 
trogram time-frequency distribution derived from the time differential impedance signal. 

Detailed Description of the Preferred Embodiments 

35 

[0026] In the various figures like reference numerals are used to refer to like elements. There is shown diagramatically 
in Fig. 1 a preferred cardiographic impedance monitoring system according to the present invention, which is indicated 
generally at 20 and is coupled to a patient P for use. System 20 preferably includes a first, outer pair of electrodes 21 
and 24, 21 being a strip electrode and 24 being a band electrode, and two pairs of parallel connected, spot-type skin 

40 electrodes, indicated generally at 22 and 23. The patient's heart is indicated diagrammatically, in phantom, at H. The 
first pair of electrodes 21 and 24 are applied to the patient's skin on opposite upper and lower sides of the heart H, 
preferably equally spaced from the heart. The pairs of parallel coupled spot electrodes 22 and 23 are applied to the 
patient on opposite upper and lower sides of the heart H respectively and between the first or outer pair of electrodes 
21 , 24. Each pair of electrodes 22 and 23 are preferably positioned on opposite lateral sides of the patient at uniform 

45 heights above and below the heart H. Again, electrode pairs 22, 23 are preferably equally spaced from or otherwise 
symmetrically positioned with respect to the heart H. A pair of conventional, electrocardiogram electrodes 25a and 25b 
are further provided. 

[0027] Preferably, all five electrodes 21-25 are coupled through a signal pickup and preprocessor 26 and then through 
an analog to digital converter 28 to a data processor 30. Preferably a color video monitor 32 and an optional hard copy 

50 printer 34 are provided under the control of the processor 30. While a separate pair of electrocardiogram electrodes 
25a and 25b is shown, it will be appreciated that the electrode pairs 22 and 23 might also be configured to generate 
an electrocardiogram signal by suitable processing of their signals. The outer pair of electrodes 21 and 24 are used 
to apply a fluctuating, preferably alternating current through the patient P between those electrodes 21 and 24 while 
electrode pairs 22 and 23 are provided to sense voltage levels on the patient P from which thoracic impedance is 

55 determined, as will be discussed. 

[0028] The general advantages tetrapolar electrode configurations, such as that of system 20, have over bipolar 
electrode configurations are that the voltage sensing electrodes 22 and 23 are substantially free of skin impedance, 
can measure impedance with less electrode-interface artifact and offer the possibility of providing a more uniform 
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current density distribution in the segment of interest in the patient. 

[0029] There is a distortion in the current density distribution in patients in the vicinity of the current electrodes 21 
and 24, known as the edge effect. In addition to edge effects, uniform current distribution is disturbed by the superpo- 
sition of a relatively highly conductive sphere, namely the heart, and the relatively low conductivity of air in the lungs 

5 in the region between electrodes 22 and 23. The current distribution has been found to be more uniform in the central 
region between the electrodes 21 and 24 the more widely spaced those electrodes are. In addition, it has been found 
that the positions of the current electrodes 21 and 24 relative to those of the voltage detecting electrodes 22 and 23 
are further important to ensure uniform current distribution. It has been found that if the distances between adjoining 
pairs of the electrodes 21 and 22, and 23 and 24 are sufficiently great, at least about 3 cm and, preferably, about 4 cm 

10 or more, the edge effect is essentially minimized and the impedance measurement becomes stable. 

[0030] The upper current electrode 21 is preferably applied to the patient's forehead while the upper voltage sensing 
current and pick-up electrodes 22 are applied to the patient's neck. This configuration has several advantages; over 
previous systems employing neck current and pickup electrodes. It is a simple configuration for the technician to re* 
member while assuring that the upper two electrodes 21 and 22 are sufficiently separate to avoid edge effects. The 

15 further spacing of the upper current electrode 21 from the heart provides a more uniformly parallel current distribution 
and equal potential lines through the patient, particularly in the area between the measuring electrodes 22 and 23. In 
addition, this arrangement minimizes the problem of attempting to locate all upper electrodes 21 and 22 on the patient's 
neck where catheters and bandages are often encountered. The forehead is usually untouched by medical devices. 
All four electrodes 21-24 are preferably mounted to the patient in parallel planes to optimize the uniformity of the current 

20 distribution and the location of the sensing electrodes 22 and 23 along the equal potential lines. 

[0031] Band electrodes have been found to provide the most parallel current lines and equal potential distributions 
through the patient and therefore should, at least theoretically, give the best results. There are several disadvantages 
to band electrodes. They may be difficult or impossible to apply to patients with extensive chest wounds because of 
the presence of dressings, tubes, etc. They require relatively long polarization time to stabilize before measurements 

25 can be taken and, therefore, may not be useable in emergency situations. They are more sensitive to motion artifact 
due to the large area of contact of the electrode. Lastly, they are generally more uncomfortable for patients when used 
for long-term monitoring. The advantages and disadvantages of spot electrodes are just the reverse of those of the 
band electrodes. 

[0032] The correlation between the results of using an all band electrode configuration and the mixed band and spot 
30 electrode configuration shown in Fig. 1 were distinctly higher (correlation coefficient equal to about 0.95) than the 
correlation between an all band and an all spot electrode configurations (correlation coefficient equal to about 0.75). 
The combined band and spot electrode configuration shown in Fig. 1 appears to have all the advantages of each type 
with almost no disadvantages of either. In particular, the mixed electrode system supplies parallel, iso-current lines 
and parallel equal potential lines generally perpendicular to the iso-current lines. The system is easy to apply to patients, 
35 even with extensive chest wounds, dressing, tubes or the like on the chest. It has relatively short polarization time and 
thus can be used in emergency situations. It is less sensitive to motion artifact due to the small area of contact provided 
by spot electrode pairs 22 and 23 which provide the impedance signal. The spot electrodes 22 and 23 are also more 
comfortable to the patient. 

[0033] The signal pick-up and preprocessor 26 is preferably a Minnesota Impedance Cardiograph (MIC) Model 304B, 
40 supplied by Sorocom, Inc., of Minnesota. The MIC 304B provides a high frequency (approximately 100 kHz), low 
amplitude (4mA RMS) alternating current at pick-ups provided for the electrodes 21 and 24. Pick-ups are also provided 
for the parallel coupled voltage sensing electrode pairs (22) and (23) and for the EKG electrodes 25a, 25b. 
[0034] The preferred MIC preprocessor generates and outputs four analog signals: the mean thoracic impedance 
signal (Z 0 ), the change in thoracic impedance signal (delta 2 or A 2), the time-derivative impedance signal (dZ/dt) and 
45 the electrocardiogram signal (EKG). The mean thoracic impedance signal, Z,,, is the impedance difference sensed 
between electrodes 22 and 23. The change in thoracic impedance signal, delta 2, is an amplification of the original 2q 
signal from which the DC component has been removed. The time-derivative impedance signal, dZ/dt, is the time 
derivative of the amplified delta Z signal. The EKG signal is conventional. The four signals are shown diagrammatically 
on monitor 32 in Fig. 1 while expanded versions of the dZ/dt and EKG signals are shown in Figs. 2A and 2B, respectively. 
so [0035] Fig. 2a depicts a portion of an exemplary EKG signal while Fig. 2b depicts a portion of an exemplary inverted 
dZ/dt signal. Fig. 2c depicts a portion of an exemplary phonocardiogram signal corresponding to the dZ/dt signal. The 
dZ/dt signal is conventionally inverted so that its maximum slopes will appear positive, thereby enabling the clinician 
to observe the cardiac event in a more familiar manner. 

[0036] Cardiac events appear in the impedance and EKG signals as perturbations or "waves". EKG waves related 
55 to each heartbeat are also noted on the EKG signal of Fig. 2a at P, Q, R, S and T. The waves which are related to a 
single: heartbeat and which appear in the dZ/dt signal are indicated in Fig. 2b at A, B, C, F, X and O. Heart sounds 
(S1 and S2) associated with the compression and relaxation, respectively, of the heart during a single beat, are overlaid 
in the dZ/dt signal of Fig. 28. In many cases, however, the identification of individual waves in either signal is not so 
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apparent. 

[0037] The four analog signals from the preprocessor 26 are passed to the analog digital converter 28. Preferably, 
the A/D converter 28 is configured for differential conversion at a sampling rate of 500 Hz for each of the four analog 
signals with a twelve bit resolution in offset binary format. The analog to digital converter may be, for example, a Data 

5 Translation Model DT 2811/PGH. 

[0038] The four digitized signals are preferably passed into an allocated memory such as a hard disk or RAMDISK 
in or associated with the processor 30 and stored in binary format for subsequent processing. The processor 30 pref- 
erably also formats the binary data signals for real time or essentially real time display on the video monitor 32. At least 
the dZ/dt and the EKG signals from preprocessor 26 are displayed for the clinician on the monitor 32 but the Zq signal 

10 and the delta Z signal are also preferably simultaneously displayed for the clinician who can check the validity of the 
various signals as they are being acquired. Preferably, system 20 processes the impedance and EKG signal data to 
determine at least the heart rate (HR) and the cardiac output (CO) of the patient in at least near real time and preferably 
displays those values as well. These and other values determined by the process may be printed out on the hard printer 
34. 

is [0039] Figs. 3a-3b depicts in block diagram form, the operational steps of the processor 30 in determining heart rate, 
cardiac output and optionally a variety of other parameters reflecting cardiac performance begins the generation of 
data. Preferably, the system 20 is configured to be operable in a predata signal processing mode in which the preproc- 
essor 26 and ADC 28 operate to pass binary data signals to the processor 30 which formats the signals for real-time 
or essentially real-time display on monitor 32. After the clinician has verified successful signal gathering and digitization 

20 through the display, the data signal processing for cardiac parameter determinations are begun. Preferably the proc- 
essor 30 will initialize the necessary calculation variables and begin storing the binary EKG and impedance signal data 
in a hard disk or RAMDISK in the binary format for subsequent processing. 

[0040] Initial data processing includes converting the stored binary impedance and EKG signal data into decimal 
integer format and reading the decimal format data into a working memory area of the processor 30 in predetermined 

25 time period blocks for processing. The EKG signal data is first processed by low-pass (LP) and high-pass (HP) filtering, 
then differentiated and thereafter non-linearly scaled or transformed. The filtered, differentiated, scaled EKG data is 
then further adaptively processed to detect the R wave peaks, which are used to identify the cardiac cycles of the 
heartbeat, and to validate the detected peaks. If valid, the processor 30 stores the time of occurrence of each validated 
peak as the peak of an R wave in the EKG signal block. 

30 [0041] In particular, preferably each consecutive five second block of the processed EKG signal data is searched 
for its maximum peak. The maximum peak is multiplied by a predetermined fractional value, preferably 0.5, to generate 
a tentative threshold. The block of data is searched to identify peaks that are greater than or equal to this threshold, 
with each consecutive pair of peaks being separated by a predetermined time interval, preferably in the range of 0.28 
to 4.5 seconds. The coefficient of variation ("C.V. n =Standard Deviation/Mean) of the peak time intervals is determined 

35 and compared to a predetermined ratio, preferably 0.3. If the determined C.V. is less than or equal to this predetermined 
ratio (0.3), the identified peaks are accepted as the R wave peaks. If the determined C.V. is larger than this ratio (0.3), 
and if the number of intervals that are less than the mean is more than the number of intervals that are greater than 
the mean, then the tentative threshold is reset by multiplying the maximum EKG value by a smaller predetermined 
fractional value, preferably 0.4. If the reverse is true, i.e., if the determined C.V. is larger than the predetermined ratio 

40 (0.3) and the number of intervals that are less than the mean is less than the number of intervals greater than the 
mean, then the tentative threshold is reset by multiplying the largest EKG value by a larger predetermined fractional 
value, preferably 0.6. The process is repeated to determine the C.V. of the newly identified intervals, as was done with 
the originally tentative threshold. If the determined C.V. still exceeds 0.3, the data block is discarded and the next block 
of processed EKG signal data is examined. 

45 [0042] After checking for control flags that may have been generated during EKG signal block processing for the 
proscribed C.V. value, the processor 30 processes the data from the impedance signals, particularly the dZ/dt signal. 
Briefly, the R wave peaks identified from the EKG signal are used to "window" or define blocks of the dZ/dt data spanning 
a single heartbeat for processing. A portion of the time-derivative impedance signal dZ/dt is selected to span one 
heartbeat of the patient. The time-derivative impedance signal preferably is smoothed by low-pass and high-pass 

50 filtering and a time-frequency distribution of the filtered data from the windowed portion of the time-derivative impedance 
signal is generated in a manner to be described. Preferably, the time-frequency distribution is thereafter analyzed, 
preferably in conjunction with the time-derivative impedance signal, to identify from the time-frequency distribution a 
time or times of occurrence of at least one, and preferably two, cardiac events in the selected portion of the time- 
derivative signal. If the processor 30 is unable to identify the waves being sought, an artifact warning control signal is 

55 generated to set a control flag. Otherwise, the pertinent parameters from the various impedance signals, as well as 
related values derived from those signals including, but not limited to, the mean impedance at the peak of the 
heartbeat contraction, dZ/dt^n, heart rate (HR), ventricular ejection time (VET), and all other parameters which may 
be identified or determined are moved into memory for storage. The process is repeated for each heartbeat identified 
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in the interval of the processed EKG signal, which was originally selected to include the several consecutive seconds 
of data with several consecutive heartbeats of the patient. When the various parameters have been determined for 
each heartbeat of the interval (QRS Empty), the processor sets control flags for any errors noted and performs any 
remaining calculations. The processor then converts the identified and determined values into floating point numbers. 
Lastly, the processor preferably displays some or all of the identified and determined values on the monitor 30 and 
prints them on the hard copier 32, if desired. 

[0043] A distinct advantage possessed by the present system over prior known systems is the filtering method em- 
ployed for processing the digitized EKG and impedance signals. 

[0044] Preferably, all filtering is digitally performed by processor 30 and the filters used are designed to employ only 
integer coefficients. This allows the filters to operate in near real-time on a relatively simple microprocessor such as 
an IBM® or IBM® compatible personal computer, which is also employed to perform the remaining processing steps. 
[0045] The all-integer filter designs of the present invention are both simple to program and fast to execute. They 
have proven to be more than capable of handling at relatively high speed the data filtering required by the present 
system. They utilize a high-level, C language without the need for assembly language programming or a co-processor. 
The filters of the invention use only a relatively small number of multiplier and additional components, all with integer 
coefficients. A brief derivation of the filters follows. 

[0046] The usual auto-regressive moving average (ARMA) system can be represented by the equation: 

y[n]=a 1 y[n-1]+...+a m y[n-m]+b 0 x[nJ+...+b k x[n-k] 
[0047] For a low-pass filter, consider a special case of the moving average (MA) system: 

y[n]=x[n]-x[n-k] 

[0048] Its transfer function is: 

H(z)=1-z" k 

where z=e sT and T is the sample interval. 
The k zeros in the z-plane are the roots of: 

1-z' k =0 

[0049] Consider the case k=12 for the z plane represented by a real abscissa (X axis) and an imaginary ordinate (Y 
axis), i.e. Z(x,y)=Z(real, imaginary). If a zero at z=(1 ,0) of the unit circle for the real and imaginary values in the z-plane, 
respectively, is cancelled, then the following low-pass filter transfer function is obtained: 

H(z)=(1-z" k )/(1-z 1 ) 
This transfer function yields the recursive relation: 

y[n]=y[n-1]+x[n]-x[n-k] 

[0050] To improve the side lobes (-14dB), second or third order zeros and poles may be taken: 

H(z)=(1-z k ) 2 /(1-z- 1 ) 2 , 

or 



H(z)=(1-z k ) 3 /(1-z 1 ) 3 



EP 0 606 301 B1 

which greatly improves the sidelobes (-27dB, -42dB). The recursive equation becomes: 

y[n]=2y[n-1]-y[n-2]+x[n]-2x[n-k]+x[n-2k]. 

5 

The total order of this system has become 2k (or 3k). For k=12, the total system needs to preserve only 24 (or 36) data 
points. 

[0051] By moving the cancelling pole z=(1,0) to z=(-1,0), a high-pass filter is obtained with the following transfer 
function: 

10 

H(z)=(1-z' k ) 2 /(1+z' 1 ) 2 

and the following recursive equation: 

15 

y[n]-2y[n-1]-y[n-2]+x[n]-2x[n-k]+x[n-2k] 

[0052] Cancelling the zeros on the unit circle at any angle q in the z-plane other than z=(1,0) and z-(-1,0) will give 
20 a band pass filter. However, if only integer coefficients are desired, some restrictions should apply since any pair of 
those zeros will result in the cosine function: 

A 2cos q ^ 1 

25 z 

Accordingly, angle q must be 60°, 90° and 120°, correspondingly located at 1/6T, 1/4T and 1/3T The bandwidth can 
then be managed by the integer k and T. 

[0053] An example of a band pass filter with a center frequency f=1/4T and nominal bandwidth NB=1/6T, would 
30 include the following transfer function and recursive equation: 

H(z)=(1-z 12 ) 2 /(1+z- 2 ) 2 

35 and 

y[n]=-2y[n-2]-y[n-4]+x[n]-2x[n-12]+x[n-24] 
40 [0054] Consider an ideal, continuous time differentiator and its respective frequency response: 

y(t) = |[x(t)] 

45 and 

H(ja)=jO 

so Since the input signal is restricted to be bandlimited, it would be satisfactory if the continuous response could be: 

H eff (jUHft, \Q\<n/T, 
55 =0, \Q\>nfT 

The corresponding discrete-time differentiator has the following frequency response and is periodic with period 2rc : 
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H(e jtt) ) =j ©rr, |to|<7t 

It can be shown that the corresponding impulse response of this discrete frequency response can be represented as: 

5 

h[n]=[njccos(nji)+sin(nK)]/n 2 jiT -»<n<oo 
which is zero for n = 0 and as follows for n not equal to zero: 

10 

h[n]=[cos(n7c)]/nT. 

[0055] For example, a 6-point differentiator would be represented as: 

15 

y[n]=(x[n+3]/3-x[n+2]/2+x[n+1]-x[n-1]+x[n-2]/2-x[n-3])n". 

[0056] A preferred, low-pass auto-regressive moving average filter is employed having a nominal bandwidth =±33. 3 
20 Hz, Sidelobes=-27dB, K=1 5, T=500 Hz/sec, Delay=1 5T, Gain=225, having the following transfer function and recursive 
formula: 

H(z)=(1-2z' 15 +z' 30 )/(1-2z' 1 +z" 2 ) 

25 

y[n]=2y[n-1]-y[n-2]+x[n]-2x[n-15]+x[n-30] 

[0057] A preferred high-pass ARM A filter is employed in the system having a bandpass =±5 Hz, and sidelobes= - 
30 1 4dB. Because of the T restriction for the specific high-pass filter, a low-pass filter is designed first and then subtracted 
by means of a linear phase delay filter, as shown in Fig. 4. The low-pass filter is represented by the following transfer 
function and recursive formula, respectively: 

3 5 HL(z)=(1-z' 100 )/(1-z' 1 ) 

y[n]=y[n-1]+x[n]-x[n-100] 

40 where K=1 00, T=500 Hz/sec, Delay=1 00T, and Gain=1 00. 

[0058] The preferred delay filter used in Fig. 4 is represented by the recursive formula: 

d[n]=100x[n-100)-y[n] 

45 

where d[n] represents the high-pass filter output samples. 

[0059] A preferred moving average (MA) differentiator used by the system 20 for differentiating the EKG signal data 
after low-pass/high-pass filtering has a linear slope bandwidth =±35 Hz and the following transfer function and recursive 
formula: 

so 

H(z)=-z^-2z" 2 +2z +2 +z +4 

55 y[n]=-x[n-4]-2x[n-2]+2x[n+2]+x[n+4] 
where K=4, T=500 Hz/sec, and Delay=4T 

[0060] Fig. 5 depicts diagramatically in block diagram form and greater detail, the sequence of steps preferably 
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followed by the processor 30 to process the digitized EKG data and to identify the R wave pulses. The EKG data stored 
in binary form is converted into digital data in blocks containing the continuous EKG signal data for a predetermined 
period of time. The blocks are selected to include several consecutive seconds of data from the EKG signal for several 
consecutive heartbeats of the patient. Five second long blocks of EKG signal data are presently preferred as such 
s blocks are sufficiently long to permit the elimination by filtering of artifacts which might be caused by arrhythmia or 
patient movement and yet are sufficiently short that they still permit the various digital filters of the system to handle 
heart rates of between about 13 and 214 beats per minute. 

[0061] Processor 30 preferably first converts the block of binary data into digital format values. It then low-pass filters 
the decimal format data, high-pass filters the low-pass filtered data, and applies the twice-filtered data to the previously 
10 described differentiator or derivative filter. The filtered, differentiated data is then non-linearly scaled, preferably by 
squaring for convenience, although raising the data by any power greater than one will provide non-linear scaling. The 
preferred low-pass ARMA filter has a nominal bandwidth equal to ±33.3 Hz to remove the components of the EKG 
signal which may be provided by external electrical interference, specifically 60 Hz AC line and other higher frequency 
interferences. The preferred high-pass ARMA filter has a band pass equal to ±5 Hz to eliminate unwanted low-frequency 
is components representing slow motion artifacts such as respiration and other physical patient movements and such 
anomalies as arrhythmia. The preferred MA differentiator has a linear low bandwidth equal to ±35 Hz to extract and 
emphasize data concerning the slope change of the EKG signal. The non-linear transformation is preferably significantly 
great to provide an R/P equal to -75 dB where R and P stand for the amplitudes of the R and P waves, respectively, 
of the EKG signal associated with each heartbeat. Squaring has proved adequate. The filtered, differentiated, non- 
20 linearly scaled five second block of EKG signal data will include a number of peak values, the greatest of which will 
correspond to the separate R wave peaks of the original EKG signal. 

[0062] According to an important aspect of the present invention, an adaptive threshold is set for each multisecond 
block of EKG data. This is done by identifying the magnitude of the maximum amplitude value of the filtered, differen- 
tiated, non-linearly scaled data of the block, setting a predetermined fraction of that magnitude, for example, .5, as a 

25 tentative threshold. All data points in the block of filtered, differentiated, non-linearly scaled data having an amplitude 
at least as great and preferably greater than that predetermined fractional threshold are identified and processed to 
determine the average of the identified amplitudes and a standard deviation of the identified amplitudes. The processor 
thereafter determines the coefficient of variation (CV) of the identified amplitudes, where CV=standard deviation/mean. 
If the coefficient of variation exceeds a predetermined value, preferably 0.3, none of the identified amplitudes are 

30 accepted or validated and that block of EKG data is preferably discarded. If the coefficient of variation of the amplitudes 
is equal to or less than the predetermined fraction, the EKG data is further processed. 

[0063] A bi-level timing signal preferably is generated from the times of occurrence of those identified amplitudes 
remaining after the foregoing processing and is output as a pulsed, heart rate signal, which can be displayed on the 
monitor 30 and/or recorded. 

35 [0064] Testing of this process on a large number of patients, including some with pacemakers and some exercising 
on a stationary bicycle, was found to provide a correct R wave detection rate for more than ninety-five percent of the 
EKG signals where a coefficient of variation of 0.3 was employed. 

[0065] Next, the processor 30 preferably processes a block of the impedance signal data overlapping the selected 
block of processed EKG data to identify cardiac events reflected in the impedance signals from which the Ventricular 
40 Ejection Times (VET) may be estimated. 

[0066] The VET plays an important role in impedance cardiography since it represents the time period between the 
opening and closing of the aortic valves during the systole-diastole cycle of the heart. It has not been clear how the 
VET should be defined with respect to the time-derivative impedance signal dZ/dt. Several alternative definitions have 
been proposed including: 

45 

(1) the distance between two zero crossing points of the dZ/dt signal; 

(2) the distance from a first zero crossing of the dZ/dt signal, before the occurrence of the "X" wave, to the "X" 
wave of the dZ/dt signal; 

(3) the distance between the 0.15* dZ/dt^ point and the X wave of the dZ/dt signal; 
so (4) the distance between the A and X waves of the dZ/dt signal; and 

(5) the distance between the B and X waves of the dZ/dt signal. 

[0067] The latter two definitions are baseline independent. The last of these definitions is presently preferred because 
the B and X waves are considered to represent the onset points for opening and closing of the valves and thus, phys- 
55 iologically correspond most closely to the end points of the theoretical interval being determined. However, this last 
definition is perhaps the definition least used because it is the most difficult one to detect. The B and X points often, if 
not usually, disappear into the noise of the basic dZ/dt signal. 

[0068] According to another important aspect of the invention, a time-frequency distribution analysis is generated to 
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extract detailed information on the transient behavior of the non-stationary dZ/dt signal. This leads to the concept of a 
mixed time-frequency signal representation. 

[0069] A general class of time-frequency distributions was introduced by Cohen and have the following form: 



where f(u) is the time signal, f*(u) is its complex conjugate, and 0 is a kernel function which represents the particular 
distribution function selected. Several distributions have been derived from this Cohen distribution including: the Rihac- 
zek distribution, the Page distribution, the Levin distribution, the Wigner distribution, the cumulative decay or attack 
spectrum and the Spectrogram distribution 

[0070] Preferably, a Spectrogram time-frequency distribution of the digitized dZ/dt signal is used by the processor 
30. The Spectrogram distribution, sometimes called the time-dependent Fast Fourier Transform (FFT), is defined gen- 
erally as: 

X[n,co] = £ f[n+mMmIc* ju)m 
m = — • 

where f[n] is the digitized dZ/dt time signal data and w[n] is a window function. Window functions particularly useful 
with Fourier transforms include rectangular, Bartlett (triangular), Hanning, Hamming, and Blackman. Preferably, a Ham- 
ming window function is used having the following equation: 

w[n] = 0.54-0.46 cos(2rcn/M), 0 < n < M 
= 0, 

otherwise. 

Variable n=0,1 M-1 and corresponds to data points of the digitized dZ/dt signal; M is the total number of data points 

(n) in the interval Rj to Rj +1 . 

[0071] When multiplied by a window function, the one-dimensional discrete signal f[n] (i.e. the digitized dZ/dt signal) 
is converted into a two-dimensional function of time variable and a frequency variable. Its time-dependent Fourier 
transform can also be interpreted as the Fourier transform of f[n+m] as viewed through the window w[m]. The window 
has a stationary origin and as n changes, the signal slides past the window so that, at each value of n, a different 
portion of the signal is viewed. 

[0072] A major advantage of time-frequency analysis utilizing a Spectrogram distribution is that for the purposes of 
the present system 20, such an analysis is possible with the computational power of the preferred microprocessor 30. 
The corresponding disadvantage of a Spectrogram analysis is that there is a trade off between time and frequency 
resolution. An underlying assumption in using a window is that the spectral characteristics of the signal being viewed 
can be considered reasonably stationary over the window's duration. The faster the spectral characteristics of the 
signal change, the shorter the window duration should be. A shorter window length provides a higher resolution in time 
changes. However, the shorter the length of the window duration, the lower the frequency resolution will be. This is 
because the resolution of narrowband components will be sacrificed as window duration decreases. To increase fre- 
quency resolution, the window length must be increased. 

[0073] The primary frequency components of the power spectrum density of typical thoracic impedance signals are 
found scattered approximately within the range of 2-40 Hz. However, the present system and method is preferably 
interested only in the contribution within the range of approximately 30-55 Hz, corresponding most closely to the fre- 
quency of opening and closing of the aortic valve. 

[0074] The Spectrogram function S is preferably calculated for each R-R interval. To be certain that the Spectrogram 
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function captures the pertinent frequency data in the dZ/dt signal associated with each heartbeat, Fast Fourier Trans- 
forms used to generate the Spectrogram are preferably determined for each heartbeat beginning at a point in time 
before the occurrence of the heartbeat. Preferably, the time period over which the Spectrogram function is calculated 
begins before the R wave peak Rj of the heartbeat in question, preferably about twenty percent of the time interval 
between the heartbeat in question (R,) and the next heartbeat(Rj +1 ). For simplicity, reference will hereafter be made 
to R-R intervals but it will be understood where it is not stated that these intervals are shifted slightly backwards in time 
for the spectrogram function computations. Preferably, the heartbeat times and intervals (R-R intervals) are defined 
by the previously derived R peaks and/or heart rate signal. The FFT is preferably calculated with the low-pass/high- 
pass filtered (smoothed) dZ/dt(n) data. 

[0075] The following variable names have been employed in Fig. 6: 

Rj : The time location of a particular R spike of the EKG signal. 

R| +1 : The time location of another R spike immediately after Rj of the same EKG signal. 

M: The number of dZ/dt sample points between the Rj and Rj +1 locations (M=Rj +1 -Rj) and is a varying number. 

D: The number of points the Rj and Rj +1 locations are shifted backwards in time. 

B: Location of B wave peak of the dZ/dt signal in the Rj to Rj +1 interval. 

X: Location of X wave peak of the dZ/dt signal in the Rj to Rj +1 interval. 

P: Location of dZ/dt^ within the Rj to Rj +1 interval. 

w[n]: A window function of data points "n". 

H[n]: Fourier transform value for a frequency value that corresponds to n. 
L: The length of an FFT process. 

m: The number of points or STEP, that a 2*L duration moves ahead (to the right) in the time domain within the 

Rj to Rj +1 cycle. Note that m«2*L and preferably equals 1. 
S: Spectrogram function. 

[0076] The signals are processed R-R cycle by R-R cycle. Again, each R-R cycle is purposely shifted backwards 
early in time by D points. D is preferably twenty percent of the total number of Rj +1 -Rj data or sample points. This is to 
assure that the A and B points of the heartbeat associated with Rj are located inside the process cycle. The preferred 
sample rate is 500 Hz/sec and per channel. 

[0077] Physiologically, within every process cycle Rj to Rj +1 , there is a peak dZ/dt value. This peak value is dZ/dt^ 
when referring to the time derivative of the impedance signal and dZ/dt^ when referring to a time derivative of an 
inverted impedance signal. In the context of the present invention, these subscripts merely identify the sign of the 
differentiated impedance signal. It is assumed that the dZ/dt max is always the positively largest one within the cycle, 
which is true physiologically. This dZ/dt^ can be easily detected from the filtered (inverted) dZ/dt signal data due to 
the fact that it has the largest positive magnitude. 

The location of dZ/dt,^ within the cycle is detected and stored as P. Variables "i" and "k" are initialized. Variable Y is 
the location pointer of the first point of each 2*L period of data copied for processing while variable "k" is the kth value 
of the copied 2*L data period. The first block of the flow chart of Fig. 6 indicates that Rj, Rj +1 , M and P values are available. 
[0078] In the second block, the first 2*L dZ/dt data points (starting from RpD) are copied into a memory buffer. Note 
that L must be an order of two and is preferably at least sixteen. In other words, it can be sixteen, thirty-two, sixty-four, 
one-hundred-twenty-eight, etc. However, the relationship 2*L<M must hold. Preferably, L equals thirty-two. 
[0079] In the third block, each of these 2*L dZ/dt data points are multiplied by a selected window function (w[n]) f 
preferably a Hamming window function, also with a length 2*L The purpose is to smooth the dZ/dt data to reduce the 
edge effect. A Rectangular window function might alternatively be used. These two window functions are described 
below: 

Rectangular 

w[n]=1 for i<n<i+2*L, 
w[n]=0, otherwise; 

Hamming 

w[n]=0.54-0.46cos(2*ji*n/M) for i<n<i+2*L, 
w[n]=0, otherwise 

[0080] In the fourth block, after the impedance data has been multiplied with the selected window function, a Fast 
Fourier Transform (FFT) is performed over these 2*L points. Note that the FFT routine requires 2*L points to return L 
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points because the result has both real and imaginary parts and each part has L points. 

[0081] Generally, the mathematical equation of the discrete Fourier transform for a frequency value that corresponds 
to n can be described as follows: 



N-l 

H[n]»21 h[p)*exp(j*2*» *n*p/N) 
p=0 

where: h[p] is the input function, N is the number of input points (i.e., N=2*L), p is the discrete variable in time domain 
(expanding from 0 to N-1), n is the discrete variable in the frequency domain and equals (0,2*L-1)) or (4r,?f) or (-L,L), 
all are equivalent. 

[0082] The result of a complete Fourier transform is a complex array: 

H[-2*L/2],...,H[-1],H[0],H[1] H[2*L72]. 

where 2*L/2 corresponds to the one-half sample frequency and 0 corresponds to zero frequency. Since the negative 
frequency part of the transform is only a mirror image of the positive part, it is usually dropped from the computation. 
[0083] Here the input function is: h[p]-w[p]*dZ/dt[p]. 

[0084] For a Hamming window function, the discrete Fourier transform is given by: 
2*L-1 

H I n] - 51 dZ/dt I P ) * ( 0 . 54 - 0 . 46cos ( 2* *cxp ( j +2*7/ *n*p/2*L) 

p-0 

where n=(-5f,^-) or (-L,L). Note that the limit of p normally might be Rj<p<Rj+2L-1. However, since every 2*L duration 
is processed individually, as far as the Fourier transform is concerned, p always begins at 0. 
[0085] The actual calculation is done by using the FFT computation algorithm to gain faster speed and the negative 
frequency part is not returned in the result. However, the principle and outcome are the same. Note that the values of 
H[n] are complex in format, i.e., real and imaginary. For -L<n<L, the set H generates two terms for every frequency, 
one for the real part, another for the imaginary part. Thus the actual result of this computation is an array of integers 
in the format: 

H[0],H[1],H[2],H[3] ( [4],...H[2*L-2],H[2*L-1]. 

where H[0] and H[1] are the real and imaginary part of the 0 frequency, H[2] and H[3] are the real and imaginary part 
of the next immediately higher frequency, H[4] and H[5] for the subsequent immediately higher frequency, and so on. 
[0086] In our application, only the magnitude of the FFT is used to generate the Spectrogram function S. That is also 
determined in the fourth block according to the equation: 



L 

FFr(n) m SQROOT (H[ 2j )*H[ 2j ] +H[ 2 j+1 ) *H [ 2j + L ] ) 

and n-0,1,...L-1. 

[0087] Note that each of the magnitude points "n" corresponds to a particular frequency and thus reflects the distri- 
bution of that frequency in this 2*L portion of the signal. The bigger the magnitude is, the larger that frequency con- 
tributes. In this case, 0-15 points are generated for each Spectrogram function point k and correspond to 0-fo/2 where 
fo is the sample frequency. Since we chose 2L-32, 0-31 represents 0-fo, and 0-15 represents The part^-fo is just 
the mirror image of 0-^-. 

[0088] In the fifth block, an individual point S[k] of the Spectrogram function is determined. A certain range of fre- 
quencies contribute to the B, C and X points on the dZ/dt signal, and that frequency range is known to be between 



13 



EP 0 606 301 B1 



about 30-55 Hz, which roughly corresponds to frequency points between n=2 and n=4 (about 31.25 and 46.875 Hz). 
This range is described in the figure as frequency range W1-W2. Similarly, points n=0 to n=2 are described in the figure 
as frequency range W0-W1 and the points n=4 to n=15 as frequency range W2-W3. 

[0089] When 0<i^P (P being the location of d^dt^), the magnitudes of points between W1-W2 are summed to- 
5 gether, and are normalized by the sum of the magnitude of points between W2-W3. When P<i<M+L, the magnitudes 
of points between W1-W2 are summed together and are normalized by the sum of the magnitude of points between 
W0-W2. S[k] contributes only one point to the time-frequency distribution at time location k. 
[0090] In the sixth block, k is increased by 1 . However, i may be increased by a step size m to the right in the time 
domain. Note that m«2*L. Thus, the second 2*L duration of the signal to be processed actually overlays with the first 
10 duration. 

[0091] Note that since every 2*L duration contributes only one point, if m, the sep size, is bigger than one, then the 
number of resultant points actually shrinks by m times in the time domain. The result can be recovered back by the 
relationship: 

15 Sfl]-S[k] for (j-k; js k+m-1 ; j++) 

Preferably, m is exactly one and no correction is needed. 

[0092] In the seventh block, k is compared with M. IF k is less than or equal to M, the whole process is repeated for 

20 the second duration for k=k+1 and i=i+1 until k>M. 

[0093] A two-dimensional representation of local frequency change amplitude in the frequency range of interest 
associated with the cardiac events in question, more approximately 30-55 Hz and, more specifically, in this example, 
31.25-46.875 Hz, which is derived from the Spectrogram distribution by integrating the Spectrogram distribution over 
that frequency range of interest or, more specifically by combining the discrete Spectrogram function contributions 

25 from the frequencies in that range of interest, is depicted in Fig. 7. Each block of the Spectrogram derived data extending 
over the time interval from Rj-.2(Rj +1 -Rj) to Rj +1 -.2(Rj +1 -Rj) is analyzed after it is generated to identify at least one wave 
in the time-derivative impedance signal data dZ/dt(n) associated with the heartbeat Rj. 

[0094] Preferably, the time of occurrence of the C wave of the heartbeat in question is first identified directly from 
the processed dZ/dt signal. The C wave has the greatest amplitude in the portion of the processed dZ/dt signal asso- 

30 ciated with the heartbeat Rj (see Fig. 2B). The Spectrogram derived data is analyzed over the time interval in question 
spanning the single heartbeat Rj to identify the maximum Spectrogram derived data amplitude value, MAX, occurring 
over that interval. Proceeding backward in time through the Spectrogram derived data from the time of occurrence of 
the C wave determined from the dz/dt signal, the first identified Spectrogram derived data value exceeding MAX/2 
encountered occurs at the occurrence of the B wave or event and is referred to as MAX1 in the program. The C wave 

35 is referred to as MAX2 in the program. The time of occurrence of the X wave is preferably determined from the interval 
between the determined B and C events. Preferably, the X wave is identified as the first peak in the Spectrogram having 
a magnitude > MAX/3 and occurring at least (B-C) seconds after the time of occurrence of the C wave peak, where B 
and C represent the times of occurrence of each of those two determined wave peaks. It is referred to as MAX3. 
[0095] Preferably, the system utilizes a VET = B - X, where B and X represent the time of occurrences of those events 

40 as determined in the above described manner the Spectrogram distribution of the dz/dt signal. Also in the preferred 
embodiment, dZ/dt^ is determined to be the absolute magnitude of the original dZ/dt signal between the times of the 
B and C events as derived from the dZ/dt and Spectrogram signals. VET and dZ/dt^ can thereafter be plugged into 
either the Xubicek or Sramek equations. 

[0096] After stroke volume (SV) is determined, for example by either the Kubicek or Sramek equations, cardiac output 
45 (CO) is determined by multiplying the stroke volume by the heart rata in beats per minute, which is preferably determined 
as described above from the EKG signal. 

[0097] In comparison parallel tests on patients determining the cardiac output values obtained by the preferred sys- 
tem, utilizing the Kubicek model with a body weight correction factor, and the invasive Swan-Ganz thermodilution 
method, correlation coefficients of about 0.8 or more were found between the cardiac outputs determined by the two 
so methods. 

[0098] Preferably, the basic Kubicek equation was modified to take into account over/underweight of the various 
patients involved by multiplying the basic Kubicek equation by a size factor F determined empirically based upon ideal 
male and female weights. 

[0099] Ideal female weight (IFW) in kg (pounds) was given by the relation: 

55 

IFW=0.534*(height in multiples of 2,54 cm (inch)) -27.36/0,454 Ideal male weight (IMW) in kg (pounds) was given 
by the relation: 

IMW=0.534*(height in multiples of 2,54 cm (inch) -1 7.36/0,454 Deviation of a patient's weight from the ideal weight 
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is given by the relation: 

deviation=(actual weight-ideal weight)/ideal weight The stroke volume calculated by the Kubicek model was 
adjusted up or down by a weight factor F, which had the following values for the following indicated deviation: 

5 

If deviation <-0.50, 

F=0. 1 0*Deviation+0.90 
If -0.5 < or equal Deviation <0.0, 

F=0.20* Deviation+1.00; 
10 If 0.0 < or equal Deviation <0.6, 

F=0.43* Deviation+1.00; 

and 

If 0.6 < or equal Deviation, 
F=0.60* Deviation+0.95. 

15 

[0100] Preferably, blood resistivity R in the Kubicek model was set equal to a constant 135 ohm-cm for male and 
112 ohm-cm for female patients. It was found that when actual value of blood resistivity was calculated by the empirical 
equation: 

20 R^53.2e {0022)HEM 

where HEM is the measured hematocrit, that significant underestimation resulted. 

[0101] Alternatively, in place of either the Kubicek or Sramek model, the processor 30 of the system 20 can be 
25 configured to implement other models including the following relation of stroke volume: 

SV=(L 2 F[weight]/C c ) (Z G /Z 2 Z 3 ] 

30 where C c is the mean thoracic conductivity, Z 2 is the impedance to ground, (absolute impedance for absolute ground) 
sensed by cne pair of voltage sensing electrodes, for example 22, and Z 3 is the absolute impedance to ground, sensed 
by the remaining pair of voltage sensing electrode, for example 23, and Zq is the mean thoracic impedance previously 
identified. Again, it may be necessary to provide a scale factor, either F(weight) previously identified or another empir- 
ically determined factor, to adjust the relation for underweight and. overweight individuals. The conductivity C c is the 

35 inverse of the resistivity. 

[0102] The above proposed relation between stroke volume and impedance is derived from taking into account the 
concentric positioning of the generally circular heart within an elliptical pair of lungs within an elliptical thoracic wall in 
the region on which the voltage sensing electrodes are placed for generation of the impedance signals. The impedance 
value Z 0 , Z 2 and Z 3 would be determined from the impedance signals at the peak of the R wave in the EKG signal, 

40 which preferably would be determined as has already been described. 

[0103] The proposed relation continues to use the mean thoracic impedance Z^ like both the Kubicek and Sramek 
models, and length L like Kubicek, but avoids the use of dZ/dt employed by both of those models. Consequently, this 
model has significantly reduced sensitivity to noise, motion artifact, respiratory cycles and ventilation effects. Again, it 
further can be adjusted to compensate for underweight and overweight individuals. 

45 [0104] While the preferred system identifies heartbeats from adaptively processed EKG signals, standard, threshold 
determined EKG heartbeat and other devices and methods for determining heartbeats can be used in impedance 
signal processing. While one preferred time-frequency distribution and one preferred window function were discussed 
in detail, other time-frequency distributions and/or window functions might be employed in the present method and 
device, particularly when greater processing capability becomes available at comparable costs. Furthermore, at least 

so some aspects of the invention lend themselves to parallel processing to improve upon the speed at which data process- 
ing and characteristic determination can be accomplished. The various aspects of the presently preferred embodiment 
are presently preferred based in large part upon data processing costs and the capability of current equipment. The 
adaptive identification of heartbeats and consequent heart rate determination by the present invention have their own 
independent value apart from their use in processing cardiac impedance signals. Other biological signals such as EEG, 

55 EMG, pulmonary pressure wave, etc. should also be adaptable to either the adaptive signal processing performed on 
the EKG signal or the time-frequency distribution processing performed on the impedance signal. 
[0105] While one preferred embodiment of the system and methods of the present invention have been disclosed 
and at least one additional alternate suggested, one of ordinary skill in the art will recognize that changes may be made 
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to the above-described embodiments of the invention without departing from the broad inventive concepts thereof. 
It is understood, therefore, that this invention is not limited to the particular embodiment disclosed but is intended to 
cover any modifications which are within the scope of the invention as defined by the appended claims. 



Claims 

1. A method of processing a time-derivative, thoracic impedance signal generated from a patient to identify events 
in the time-derivative impedance signal associated with beats of the patient's heart, characterized by the steps of: 

10 

- generating a time-frequency distribution of data from the time-derivative impedance signal for a selected por- 
tion of the time-derivative impedance signal spanning a single heartbeat of the patient; and 

identifying a time of occurrence of at least one cardiac event in the selected portion of said time-derivative 
15 impedance signal from the time-frequency distribution. 

2. The method of claim 1 , characterized in that the identifying step comprises the steps of: 

determining a maximum amplitude value of the time-frequency distribution; and 

20 

identifying an amplitude value of the time-frequency distribution at least as great as a predetermined fractional 
value of the maximum amplitude value. 

3. The method of claim 2, characterized in that the amplitude value identifying step comprises the steps of: 

25 

initially identifying a time of occurrence of a maximum amplitude value in the selected portion of the time- 
derivative impedance signal; and 

identifying the first amplitude value of the time-frequency distribution at least as great as the predetermined 
30 fractional value occurring in the time-frequency distribution most immediately preceding the time of the occur- 

rence of the maximum amplitude value of the time-derivative impedance signal. 

4. The method of claim 3 further characterized by the steps of determining stroke volume of the patient during the 
single heartbeat from the difference between the time of occurrence of the first amplitude value of the time-fre- 

35 quency distribution and the time of occurrence of the maximum amplitude value of the time-derivative signal. 

5. The method of claim 3 further characterized by the step of identifying a second amplitude value of the time- 
frequency distribution at least as great as the predetermined fractional value and following the time of occurrence 
of the maximum amplitude value of the time-derivative signal by at least a predetermined time period. 

40 

6. The method of claim 5 further characterized by the step of determining stroke volume of the patient during the 
single heartbeat from the difference between the time of occurrence of the first amplitude value of the time-fre- 
quency distribution and the time of occurrence of the second amplitude value of the time-frequency distribution. 

45 7. The method of claim 2, characterized in that the amplitude value identifying step comprises the steps of: 

initially identifying a time of occurrence of a maximum amplitude value in the selected portion of the time- 
derivative impedance signal; and 

50 - identifying a second amplitude value of the time-frequency distribution at least as great as the predetermined 

fractional value and following the time of occurrence of the maximum amplitude value of the time-derivative 
signal by at least a predetermined time period. 

8. The method of claim 1 , characterized in that the steps of generating the time-frequency distribution comprises 
55 generating a Cohen derived time-frequency distribution of the selected portion of the time-derivative signal. 

9. The method of claim 8, characterized in that the step of generating the time-frequency distribution further com- 
prises generating the distribution with a window function. 
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10. The method of claim 1 further characterized by the initial steps of: 

differentiating an electrocardiogram signal from the patient with respect to time; 

5 - scaling the differentiated electrocardiogram signal non-li nearly to emphasize amplitude peaks in the differen- 

tiated electrocardiogram signal; 

identifying a maximum amplitude value of the scaled signal occurring in an interval including several consec- 
utive seconds of the scaled electrocardiogram signal; 

10 

identifying each amplitude value of the scaled signal in the selected interval at least as great as a predetermined 
fraction of the maximum amplitude value; and 

generating the time-frequency distribution from a portion of the time-derivative signal spanning the time of 
15 occurrence of only one of the identified amplitude values of the scaled signal. 

11. The method of claim 1 further characterized by the step of initially filtering the time-derivative signal to remove 
low frequency components from patient movements and higher frequency components from sources of electrical 
interference external to the patient and wherein the generating step comprises generating the time-frequency 

20 distribution from the filtered time-derivative signal. 

12. A method of claim 1 1 , characterized in that the filtering step comprises removing frequency components of about 
5 hertz or less and about 50 hertz or more from the time-derivative impedance signal. 

25 13. The method of claim 1 characterized by the steps of: 

differentiating an electrocardiogram signal data with respect to time; 

scaling the differentiated signal data nonlineariy to emphasize amplitude peaks in the differentiated signal data; 

30 

identifying a maximum amplitude value of the scaled signal data derived from an interval of the electrocardi- 
ogram signal selected to include several consecutive seconds of data for several consecutive heartbeats of 
the patient; 

35 identifying each amplitude value of the scaled signal data in the selected interval at least as great as a pre- 

determined fraction of the maximum amplitude value; and 

generating a coefficient of variation for the amplitude values of the scaled signal identified from the selected 
interval in the previous step. 

40 

14. The method of claim 13 further characterized by the step of transmitting a bi-level timing signal spanning the 
selected interval of the scaled signal and including level transitions corresponding in time to the occurrence of 
each of the identified amplitude values of the scaled signal during the selected interval. 

45 15. The method of claim 13 further characterized by the step of filtering the electrocardiogram signal to remove low 
frequency components of the signal from patient movement and high frequency components of the signal from 
sources of external electrical interference before the differentiating step. 

16. The method of claim 1 5, characterized in that the step of filtering comprises passing data of the electrocardiogram 
50 signal through a high-pass filter to eliminate frequency components of about five hertz or less from the signal data, 

17. The method of claim 1 6, characterized in that the step of filtering further comprises passing the high-pass filtered 
data through a low-pass filter to at least eliminate components of about fifty hertz or more from the high-pass 
filtered data. 

55 

18. The method of claim 16, characterized in that the step of passing the high-pass filtered data further comprises 
passing said filtered data through a low-pass filter to at least eliminate components of about forty hertz or more 
from the high-pass filtered data. 
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19. The method of claim 15, characterized in that the electrocardiogram signal data is digitized and the filtering step 
comprises passing the digitized electrocardiogram data through an all-integer coefficient filter. 

20. The method of claim 13 further characterized by the steps of: 

5 

comparing the coefficient of variation to a predetermined valued; and 

transmitting a bi-level timing signal spanning the selected interval and including level transitions corresponding 
in time to the time of occurrence of each of the identified amplitude values of the scaled signal during the 
10 selected interval only if the coefficient of the variation for the interval is at least less than a predetermined 

fractional value. 

21. Apparatus for implementing the method of claim 1 or one of the claims depending from claim 1 , characterized by 

15 - the provision of a first pair of electrodes (21 , 24) on opposing upper and lower sides of the patient's heart (H) 

to which a fluctuating excitation current is applied; 

- the provision of a second pair of electrodes (22, 23) on opposing sides of the patient's heart (H) between the 
first pair of electrodes (21 , 24) in order to generate a first signal related to impedance of the patient at a first 
20 electrode of the second pair and a second signal related to impedance of the patient of a second electrode of 

the second pair; 

means (26, 28, 30) for identifying an impedance from each of the first and second signals occurring simulta- 
neously with one heartbeat of the patient including means for determining the time-derivative of said signals; 
25 and 

means (30) for multiplying together the identified impedances of each of the first and second signals to estimate 
a cardiac parameter of the patient's heart. 

30 

Patentanspruche 

1. Verfahren zur Verarbeitung eines zeitlich abgeleiteten, von einem Patienten erzeugten Thoraximpedanzsignals 
zum Identifizieren von Ereignissen in dem zeitlich abgeleiteten Impedanzsignal, die Schlagen des Herzen des 

35 Patienten zugeordnet sind, gekennzeichnet durch die folgenden Schritte: 

Erzeugen einer Zeit-Frequenz-Verteilung von Daten aus dem zeitlich abgeleiteten Impedanzsignal fur einen 
ausgewahlten Teil des zeitlich abgeleiteten Impedanzsignals, der einen einzelnen Herzschlag des Patienten 
uberspannt; und 

40 

Identifizieren eines Zeitpunkts des Auftretens mindestens eines Herzereignisses im ausgewahlten Teil des 
zeitlich abgeleiteten Impedanzsignals aus der Zeit-Frequenz-Verteilung. 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daft der Identifizierungsschritt folgende Schritte umfalit: 

45 

Bestimmen eines gro&ten Amplitudenwerts der Zeit-Frequenz-Verteilung; und 

Identifizieren eines Amplitudenwerts der Zeit-Frequenz-Verteilung, der mindestens so groft ist wie ein vorbe- 
stimmter Bruchwert des groftten Amplitudenwerts. 

50 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daft der Amplitudenwertidentifizierungsschritt die folgen- 
den Schritte umfadt: 

anfangliches Identifizieren eines Zeitpunkts des Auftretens eines groBten Amplitudenwerts in dem ausgewahl- 
55 ten Teil des zeitlich abgeleiteten Impedanzsignals; und 

Identifizieren des ersten Amplitudenwerts der Zeit-Frequenz-Verteilung, der mindestens so grofi ist wie der 
vorbestimmte Bruchwert, der in der Zeit-Frequenz-Verteilung dem Zeitpunkt des Auftretens des groftten Am- 
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plitudenwerts des zeitlich abgeleiteten Impedanzsignals unmittelbar vorausgeht. 

4. Verfahren nach Anspruch 3, weiterhin gekennzeichnet durch den Schritt des Bestimmens des Schlagvolumens 
des Patienten wahrend des einzelnen Herzschlags aus der Differenz zwischen dem Zeitpunkt des Auftretens des 

5 ersten Amplitudenwerts der Zeit-Frequenz-Verteilung und dem Zeitpunkt des Auftretens des groftten Amplituden- 

werts des zeitlich abgeleiteten Signals. 

5. Verfahren nach Anspruch 3, weiterhin gekennzeichnet durch den Schritt des Identifizierens eines zweiten Am- 
plitudenwerts der Zeit-Frequenz-Verteilung, der mindestens so groB ist wie der vorbestimmte Bruchwert und urn 

10 mindestens eine vorbestimmte Zeitperiode dem Zeitpunkt des Auftretens des groBten Amplitudenwerts des zeitlich 

abgeleiteten Signals folgt. 

6. Verfahren nach Anspruch 5, weiterhin gekennzeichnet durch den Schritt des Bestimmens des Schlagvolumens 
des Patienten wahrend des einzelnen Herzschlags aus der Differenz zwischen dem Zeitpunkt des Auftretens des 

15 ersten Amplitudenwerts der Zeit-Frequenz-Verteilung und dem Zeitpunkt des Auftretens des zweiten Amplituden- 

werts der Zeit-Frequenz-Verteilung. 

7. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB der Amplitudenwertidentifizierungsschritt die folgen- 
den Schritte umfaBt: 

20 

anfang Itches Identifizieren eines Zeitpunkts des Auftretens eines groBten Amplitudenwerts in dem ausgewahl- 
ten Teil des zeitlich abgeleiteten Impedanzsignals; und 

Identifizieren eines zweiten Amplitudenwerts der Zeit-Frequenz-Verteilung, der mindestens so groB ist wie 
25 der vorbestimmte Bruchwert und um mindestens eine vorbestimmte Zeitperiode dem Zeitpunkt des Auftretens 

des groBten Amplitudenwerts des zeitlich abgeleiteten Signals folgt. 

8. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daft der Schritt des Erzeugens der Zeit-Frequenz-Vertei- 
lung das Erzeugen einer nach Cohen abgeleiteten Zeit-Frequenz-Verteilung des ausgewahlten Teils des zeitlich 

30 abgeleiteten Signals umfaBt. 

9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, daB der Schritt des Erzeugens der Zeit-Frequenz-Vertei- 
lung weiterhin das Erzeugen der Verteilung mit einer Fensterfunktion umfafit. 

35 10. Verfahren nach Anspruch 1 , weiterhin gekennzeichnet durch die folgenden anfanglichen Schritte: 

Differenzieren eines Elektrokardiogrammsignals von dem Patienten nach der Zeit; 

nichtlineares Skalieren des differenzierten Elektrokardiogrammsignals, um Amplitudenspitzen in dem diffe- 
40 renzierten Elektrokardiogrammsignal hervorzuheben; 

Identifizieren eines groBten Amplitudenwerts des skalierten Signals, der in einem Zeitraum auftritt, der mehrere 
aufeinanderfolgende Sekunden des skalierten Elektrokardiogrammsignals enthalt; 

45 - Identifizieren jedes Amplitudenwerts des skalierten Signals in dem ausgewahlten Zeitraum, der mindestens 

so groB. ist wie ein vorbestimmter Bruchteil des groBten Amplitudenwerts; und 

Erzeugen der Zeit-Frequenz-Verteilung aus einem Teil des zeitlich abgeleiteten Signals, der den Zeitpunkt 
des Auftretens von nur einem der identifizierten Amplitudenwerte des skalierten Signals uberspannt. 

50 

11. Verfahren nach Anspruch 1, weiterhin gekennzeichnet durch den Schritt des anfanglichen Fi Items des zeitlich 
abgeleiteten Signals zum Entfernen von niederfrequenten Anteilen von Patientenbewegungen und hoherfrequen- 
ten Anteilen von Quellen elektrischer Interferenzen auBerhalb des Patienten und wobei der Erzeugungsschritt das 
Erzeugen der Zeit-Frequenz-Verteilung aus dem gefilterten zeitlich abgeleiteten Signal umfafit. 

55 

12. Verfahren nach Anspruch 11, dadurch gekennzeichnet, daB der Filterungsschritt das Entfernen von Frequenz- 
anteilen von etwa 5 Hz Oder weniger und etwa 50 Hz Oder daruber aus dem zeitlich abgeleiteten Impedanzsignal 
umfaBt. 
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13. Verfahren nach Anspruch 1, gekennzeichnet durch die folgenden Schritte: 

Differenzieren von Elektrokardiogrammsignaldaten nach der Zeit; 

5 nichtlineares Skalieren der differenzierten Signaldaten, urn Amplitudenspitzen in den differenzierten Signal- 

daten hervorzuheben; 

Identifizieren eines groBten Amplitudenwerts der skalierten Signaldaten, die von einem Interval I des Elektro- 
kardiogrammsignals abgeleitet sind, das so ausgewahlt ist, dali es mehrere aufeinanderfolgende Sekunden 
10 von Daten fur mehrere aufeinanderfolgende Herzschlage des Patienten enthalt; 

Identifizieren jedes Amplitudenwerts der skalierten Signaldaten in dem ausgewahlten Interval), der mindestens 
so groR ist wie ein vorbestimmter Bruchteil des grofcten Amplitudenwerts; und 

15 Erzeugen eines Variationskoeffizienten fur die aus dem ausgewahlten intervall in dem vorausgegangenen 

Schritt identifizierten Amplitudenwerte des skalierten Signals. 

14. Verfahren nach Anspruch 1 3, weiterhin gekennzeichnet durch den Schritt des Ubertragens eines zweipegeligen 
Zeitsteuerungssignals, das das ausgewahlte Intervall des skalierten Signals uberspannt und Pegelubergange ent- 

20 halt, die zeitlich dem Auftreten jedes der identifizierten Amplitudenwerte des skalierten Signals wahrend des aus- 

gewahlten Interval Is entsprechen. 

15. Verfahren nach Anspruch 13, weiterhin gekennzeichnet durch den Schritt des Filtems des Elektrokardiogramm- 
signals zum Entfernen von niederfrequenten Anteilen des Signals von Patientenbewegungen und hochfrequenten 

25 Anteilen des Signals von Quellen au&erer elektrischer Interferenzen vor dem Differenzierungsschritt. 

16. Verfahren nach Anspruch 15, dadurch gekennzeichnet, dad der Schritt des Filterns umfalit, Daten des Elektro- 
kardiogrammsignals durch ein Hochpaftfilter zu schicken, urn aus den Signaldaten Frequenzanteile von etwa 5 
Hz Oder weniger zu entfernen. 

30 

17. Verfahren nach Anspruch 16, dadurch gekennzeichnet, daft der Schritt des Filterns weiterhin umfa&t, die hoch- 
pa&gefilterten Daten durch ein Hefpaftfilterzu schicken, urn aus den hochpaligefilterten Daten mindestens Anteile 
von etwa 50 Hz Oder dariiber zu entfernen. 

35 18. Verfahren nach Anspruch 16, dadurch gekennzeichnet, daft der Schritt des Filtems der hoch pa fig efitterten Daten 
weiterhin umfaBt, die gefilterten Daten durch ein Tlefpalifilter zu schicken, urn aus den hochpaligefilterten Daten 
mindestens Anteile von etwa 40 Hz oder dariiber zu entfernen. 

19. Verfahren nach Anspruch 15, dadurch gekennzeichnet, daft die Elektrokardiogrammsignaldaten digitalisiert wer- 
40 den und der Filterungsschritt umfaftt, die digitalisierten Elektrokardiogrammdaten durch ein Ganzzahl-Koeffizien- 

ten-Filter zu schicken. 

20. Verfahren nach Anspruch 13, weiterhin gekennzeichnet durch die folgenden Schritte: 

45 - Vergleichen des Variationskoeffizienten mit einem vorbestimmten Wert; und 

Ubertragen eines zweipegeligen Zeitsteuerungssignals, das das ausgewahlte Intervall uberspannt und Pege- 
lubergange enthalt, die zeittich dem Zeitpunkt des Auftretens jedes der identifizierten Amplitudenwerte des 
skalierten Signals wahrend des ausgewahlten Intervalls entsprechen, nurdann, wenn der Variationskoeffizient 
so fur das Intervall mindestens unter einem vorbestimmten Bruchwert liegt. 

21. Vorrichtung zum Impiementieren des Verfahrens von Anspruch 1 Oder einem der von Anspruch 1 abhangigen 
Anspruche, gekennzeichnet durch: 

55 - die Bereitstellung eines ersten Paars von Elektroden (21 , 24) an einer oberen und entgegengesetzten unteren 

Seite des Herzens (H) des Patienten, an die ein schwankender Erregerstrom angelegt wird; 

die Bereitstellung eines zweiten Paars von Elektroden (22, 23) an entgegengesetzten Seiten des Herzens (H) 
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des Patienten zwischen dem ersten Paar von Elektroden (21 , 24), um ein erstes Signal, das zu der Impedanz 
des Patienten an einer ersten Elektrode des zweiten Paars in Beziehung steht, und ein zweites Signal, das 
zu der Impedanz des Patienten an einer zweiten Elektrode des zweiten Paars in Beziehung steht, zu erzeugen; 

- ein Mittel (26, 28, 30) zum Identifizieren einer Impedanz aus jedem des ersten und zweiten Signals, die gleich- 
zeitig mit einem Herzschlag des Patienten auftreten, einschlielilich einem Mittel zum Bestimmen der zeitlichen 
Ableitung der Signale; und 

ein Mittel (30) zum Zusammenmultiplizieren der identifizierten Impedanzen jedes des ersten und zweiten Si- 
gnals zum Abschatzen eines Herzparameters des Herzen des Patienten. 



Revendications 

1. Methode de traitement d'un signal d'impedance de thoracle variant dans le temps genere depuis un patient pour 
identifier les evenements dans le signal d'impedance variant dans le temps associes avec les battements du coeur 
du patient, caract6risee par les etapes consistant a : 

generer une distribution temps-frequence des donnees du signal d'impedance variant dans le temps pendant 
une portion choisie du signal d'impedance variant dans le temps couvrant un seul battement du coeur du 
patient ; et 

identifier un moment d'apparition d'au moins un evenement cardiaque dans la portion choisie dudit signal 
d'impedance variant dans le temps a partir de la distribution temps-frequence. 

2. M6thode selon la revendication 1 , caracterisee en ce que I'etape d'identification comprend les etapes consistant 
a: 

determiner une valeur d'amplitude maximale de la distribution temps-frequence ; et 

identifier une valeur d'amplitude de la distribution temps-frequence au moins aussi 6lev6e qu'une valeur frac- 

tionnelle predeterminee de la valeur d'amplitude maximale. 

3. Methode selon la revendication 2, caracterisee en ce que I'etape d'identification de la valeur de ('amplitude com- 
prend les etapes consistant a : 

identifier initialement un moment d'apparition d'une valeur d'amplitude maximale dans la portion choisie du 
signal d'impedance variant dans le temps ; et 

identifier la premiere valeur d'amplitude de la distribution temps-frequence au moins aussi elevee que la valeur 
fractionnelle predeterminee qui apparait dans la distribution temps-frequence et qui precede le plus directe- 
ment le moment d'apparition de la valeur maximale d'amplitude du signal d'impedance variant dans le temps. 

4. Methode selon la revendication 3, egalement caracterisee par les etapes consistant a determiner le deplacement 
volumique du patient pendant I'unique battement du coeur a partir de la difference entre le moment d'apparition 
de la premiere valeur d'amplitude de la distribution temps-frequence et le moment d'apparition de la valeur d'am- 
plitude maximale du signal variant dans le temps. 

5. Methode selon la revendication 3, egalement caract6risee par I'etape consistant a identifier une deuxieme valeur 
d'amplitude de la distribution temps-frequence au moins aussi elevee que la valeur fractionnelle predeterminee 
et suivant le moment d'apparition de la valeur d'amplitude maximale du signal variant dans le temps d'au moins 
une periode de temps predeterminee. 

6. Methode selon la revendication 5, egalement caracterisee par I'etape consistant a determiner le deplacement 
volumique du patient pendant I'unique battement de coeur a partir de la difference entre le moment d'apparition 
de la premiere valeur d'amplitude de la distribution temps-frequence et le moment d'apparition de la deuxieme 
valeur d'amplitude de la distribution temps-frequence. 

7. Methode selon la revendication 2, caracterisee en ce que I'etape d'identification de la valeur de ('amplitude com- 
prend les etapes consistant a : 
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identifier initialement un moment d'apparition d'une valeur d'amplitude maximale dans la portion choisie du 
signal d'impedance variant dans le temps ; et 
- identifier une deuxfeme valeur d'amplitude de la distribution temps-frequence au moins aussi elevee que la 
valeur fractionnelle pred&erminee et suivant le moment d'apparition de la valeur d'amplitude maximale du 
5 signal variant dans le temps d'au moins une periode de temps predetermine. 

8. M6thode selon la revendication 2, caract6ris6e en ce que les Stapes de generation de la distribution temps- 
frequence comprennent la generation d'une distribution temps-frequence de Cohen deriv6e de la portion choisie 
du signal variant dans le temps. 

10 

9. M6thode selon la revendication 8, caracttrisie en ce que I'etape de generation de la distribution temps-frequence 
comprend egalement la generation de la distribution avec une fonction de fenetre. 

10. Methode selon la revendication 1 , egalement caract6ris6e par les eta pes initiales consistant a : 

15 

differencier un signal d'eiectrocardiogramme du patient en fonction du temps ; 

demultiplier la non-linearite du signal d'eiectrocardiogramme differencie pour accentuer les cretes d'amplitude 
dans le signal d'eiectrocardiogramme differencie ; 

identifier une valeur d'amplitude maximale du signal demultiplie se produisant dans un intervalle comprenant 
20 plusieurs secondes consecutives du signal d'eiectrocardiogramme demultiplie ; 

identifier chaque valeur d'amplitude du signal demultiplie dans I'intervalle seiectionne au moins aussi elevee 
qu'une fraction pr6determin6e de la valeur maximale de I'amplitude ; et 

gen6rer la distribution temps-frequence a partir d'une portion du signal variant dans le temps couvrant le 
moment d'apparition d'au moins une des valeurs d'amplitude identifiee du signal demultiplie. 

25 

11. Methode selon la revendication 1, egalement caract6risee par I'etape consistant a filtrer initialement le signal 
variant dans le temps afin de supprimer les composantes a basse frequence issues des mouvements du patient 
et les composantes a frequence plus elevee issues des sources ^interferences eiectriques externes au patient 
et dans laquelle I'etape de generation comprend la generation de la distribution temps-frequence a partir du signal 

30 variant dans le temps filtre. 

12. Methode selon la revendication 11, caract6ris6e en ce que I'etape de filtrage comprend la suppression des com- 
posantes dont la frequence est d'environ 5 Hertz ou moins et d'environ 50 Hertz ou plus du signal d'impedance 
variant dans le temps. 

35 

13. Methode seion la revendication 1, caract6ris6e par les etapes consistant a : 

differencier les donnees d'un signal d'eiectrocardiogramme en fonction du temps ; 
demultiplier la non-linearite des donnees du signal differencie afin d'accentuer les cretes d'amplitude dans les 
40 donnees du signal differencie ; 

identifier une valeur d'amplitude maximale des donnees du signal demultiplie derivee d'un intervalle du signal 
d'eiectrocardiogramme choisi de maniere a inclure plusieurs secondes consecutives de donnees pour plu- 
sieurs battements cardiaques consecutifs du patient ; 

identifier chaque valeur d'amplitude des donnees du signal demultiplie dans Tintervalle seiectionne au moins 
45 aussi elevee qu'une fraction predeterminee de la valeur d'amplitude maximale; et 

g6n6rer un coefficient de variation des valeurs d'amplitude du signal demultiplie identifiees a partir de I'inter- 
valle seiectionne dans I'etape precedente. 

14. Methode selon la revendication 13, egalement caract6risee par I'etape consistant a transmettre un signal de 
so temporisation a deux niveaux couvrant I'intervalle seiectionne du signal demultiplie et comprenant des transitions 

de niveau correspondant au moment d'apparition de chacune des valeurs d'amplitude identifiees du signal demul- 
tiplie pendant I'intervalle seiectionne. 

15. Methode selon la revendication 13, egalement caract6ris6e par I'etape consistant a filtrer le signal d'eiectrocar- 
55 diogramme pour eiiminer les Composantes a basse frequence du signal issues des mouvements du patient et les 

composantes a frequence plus elevee issues des sources ^interferences eiectriques externes avant I'etape de la 
differenciation. 
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16. Methods selon la revendication 15, caract6risee en ce que I'etape consistant a filtrer com pre nd le transfert des 
donnees du signal d'electrocardiogramme au travers d'un filtre passe-haut pour eliminer des donnees du signal 
les composantes dont la frequence est d'environ cinq Hertz ou moins. 

17. Methode selon la revendication 16, caracterisee en ce que I'etape consistant a filtrer comprend egalement le 
transfert des donnees filtrees par un filtre passe-haut au travers d'un filtre passe-bas pour au moins eliminer des 
donnees filtrees par un filtre passe-haut les composantes dont la frequence est d'environ cinquante Hertz ou plus. 

18. Methode selon la revendication 16, caracterisee en ce que I'etape de transfert des donnees filtrees par un filtre 
passe-haut comprend egalement le transfert desdites donnees filtrees au travers d'un filtre passe-bas pour au 
moins eliminer des donnees filtrees par un filtre passe-haut les composantes dont la frequence est d'environ 
quarante Hertz ou plus. 

19. Methode selon la revendication 15, caracterisee en ce que les donnees du signal d'electrocardiogramme sont 
numerisees et I'etape de filtrage comprend le transfert des donnees d'electrocardiogramme numerisees au travers 
d'un filtre a coefficient entier. 

20. Methode selon la revendication 13, egalement caracterisee par les etapes consistant a : 

- comparer le coefficient de variation avec une valeur predeterminee ; et 

- transmettre un signal de temporisation a deux niveaux couvrant I'intervalle selectionne et comprenant les 
transitions de niveau correspondant dans le temps au moment d'apparition de chacune des valeurs d'amplitude 
identifies du signal demultiplie pendant I'intervalle selectionne seulement si le coefficient de variation de 
I'intervalle est au moins inferieur a une valeur fractionnelle predeterminee. 

21 . Appareil pour mettre en oeuvre la methode de la revendication 1 ou Tune des revendications dependantes de la 
revendication 1, caracterise par 

la presence d'une premiere paire d'electrodes (21 , 24) sur les cotes superieur et inferieur opposes du coeur 
du patient (H) auxquelles est applique un courant d'excitation fluctuant ; 

la presence d'une deuxieme paire d'electrodes (22, 23) sur des cotes opposes du coeur du patient (H) entre 
la premiere paire d'electrodes (21, 24) afin de generer un premier signal en rapport avec I'impedance du 
patient sur une premiere electrode de la deuxieme paire et un deuxieme signal en rapport avec ('impedance 
du patient sur une deuxieme electrode de la deuxieme paire ; 

- des moyens (26, 28, 30) pour identifier une impedance a partir de chacun des premier et deuxieme signaux 
apparaissant simultanement avec un battement cardiaque du patient comprenant des moyens pour determiner 
la variation dans le temps desdits signaux ; et 

- des moyens (30) pour multiplier ensemble les impedances identifies de chacun des premier et deuxieme 
signaux afin d'estimer un parametre cardiaque du coeur du patient. 
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